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Abstract 

KgKg production in two-photon collisions has been studied using a 397.6 ftr 1 data 
sample collected with the Belle detector at the KEKB e + e~ collider. For the first 
time the cross sections are measured in the two-photon center-of-mass energy range 
between 2.4 GeV and 4.0 GeV and angular range |cos#*| < 0.6. Combining the 
results with measurements of 77 — > K + K~ from Belle, we observe that the cross 
section ratio a (Kg Kg) /a(K + K~) decreases from ~0.13 to ~0.01 with increasing 
energy. Signals for the XcO and Xc2 charmonium states are also observed. 

Key words: Two-photon collisions; Mesons; QCD; Charmonium 
PACS: 12.38.Qk;13.85.Lg;13.66.Bc;13.25.Gv 
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1 Introduction 



Exclusive processes with hadronic final states in two-photon collision are an 
excellent probe to test various model calculations motivated by perturbative 
and non-perturbative QCD. As shown by Brodsky and Lepage (BL) [1], at 
sufficiently large two-photon center-of-mass energy and momentum trans- 
fer from the initial photon to the produced meson t, the leading term of the 
amplitude for the process 77 — > MM, where M denotes a meson, can be 
expressed as a hard scattering amplitude for 77 — > qqqq times the leading 
term meson electromagnetic form factor. For mesons with zero helicity their 
calculation gives the following dependence on s and scattering angle 9*: 



J --lead I plead / \ 1 2 

16 ^ A/ {S)l x 



d\cos9*\ s 

[( ei -e 2 ) 2 ] 2 , 2(e 1 e 2 )[(e 1 -e 2 ) 2 ; 



s(r) + 2( ei e 2 )y(0*)}, (i) 



cos 2 9*) 2 1 — cos 2 9* 

where e\ and e 2 are the quark charges (i.e., mesons have charges ±(ei — e 2 )), 
and explicit forms of the leading term meson form factor F^ ad (s) (Fj^ ad (s) ~ 
1/s at s — > oo) and the function g(9*) can be found in Refs. [1,2]. Eq. (1) 
implies that the angular distribution of neutral meson pairs, unlike that for 
charged meson pairs which is dominated by ~ sin -4 9* terms, is directly de- 
termined by the shape of g{9*) and the value of F^ d (s). Later, Benayoun and 
Chernyak (BC) [2] used a factorization hypothesis similar to the BL calcu- 
lation but further improved the treatment of the effects of SU(3) symmetry 
breaking; their predictions appeared to be in good agreement with the subse- 
quent measurements of 77 — > tx + ti~ and 77 — > K + K~ [3,4]. 

Recently, Diehl, Kroll and Vogt (DKV) [5] considered the consequences of 
the assumption that at intermediate energies the amplitudes for the process 
77 — > MM are dominated by so-called handbag contributions. The handbag 
amplitude is expressed as the product of an amplitude for the hard 77 — * qq 
subprocess times an unknown form factor RmJj{ s ) describing the soft transi- 
tion from the qq to the meson pair. In [5] the differential cross section is given 
by 

do , „ 87ra 2 1 ,„ . .„ 

"(77 - MM) = — \R MjI (s)\\ (2) 



d\cos9*\ s sin 9* 

where the meson annihilation form factor R m jj(s) is not calculated in Ref. [5] 
but is instead obtained by fitting the data; the magnitude of R M j^(s) for 
different mesons can be linked by using SU(3) and isospin symmetry. The 
validity of this approach has recently been criticized in Ref. [6]. 
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Earlier, the Belle Collaboration performed a high-statistics measurement of 
the cross sections for the processes 77 — > and 77 — > K + K~ [4] in the 

W(=yfs) range 2.4 GeV < W < 4.1 GeV. Analysis of the data showed 
that in this W range the W dependence of the cross section is consistent with 
that predicted by the leading term QCD calculations [1,2]. Here we report a 
measurement of the cross section for 77 — > K° S K° S at 2.4 GeV < W < 4.0 GeV 
and I cos 6* | < 0.6 with a data sample of 397.6 fb _1 collected at or near 
the T(4S f ) resonance, accumulated with the Belle detector [7] at the KEKB 
asymmetric-energy e + e~ collider [8]. This measurement can provide important 
information that complements previous studies and sheds light on how the 
two-photon mass and angular distributions of such cross sections depend on 
the flavor of the produced mesons. 

The Belle detector is a large-solid-angle magnetic spectrometer. Momenta of 
charged tracks are measured with a central drift chamber (CDC), located in a 
uniform 1.5 T magnetic field which surrounds the interaction point (IP) and 
subtends the polar angle range 17° < #i a b < 150°, where #i a b is a scatter- 
ing angle in the laboratory frame. The trajectories of the charged tracks near 
the interaction point are provided by the CDC and the silicon vertex detec- 
tor (SVD). Energy measurement of electromagnetically interacting particles 
is performed in an electromagnetic calorimeter (ECL) made up of CsI(Tl) 
crystals. The detector is described in detail elsewhere [7]. 



2 Event Selection 

Exclusive KgKg pairs are produced in quasi-real two-photon collisions through 
the process e + e~ — > e + e~77 — > e + e~K s K s , where the scattered e + and e~ 
are lost down the beampipe, and only the two Kg mesons are detected. 

We select 77 — > KgKg candidate events in two stages. At stage I the following 
requirements are applied: 

• exactly four charged tracks with zero net charge of which at least two have 
Pt > 0.3 GeV/c, dr < 1 cm, \dz\ < 5 cm, where pt is the transverse mo- 
mentum in the laboratory frame and dr and dz are the radial and axial 
coordinates of the point of closest approach of the track to the nominal IP, 
respectively, and the z-axis is the direction opposite to the positron beam 
axis; 

• the sum of the magnitudes of the momenta of all tracks, Hp, and the total 
energy deposit in the ECL are less than 6 GeV/c and 6 GeV, respectively; 

• the invariant mass of these four tracks is less than 4.5 GeV/c 2 , and the 
missing mass squared of the event is greater than 2 GeV 2 /c 4 ; 
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At stage II pairs of oppositely charged tracks without particle identification 
are used to reconstruct Kg — > ix + ti~ decays. To distinguish 77 — > KgKg 
events from other four-track background sources such as 77 — > 2(tt + tt~), 77 — > 
2(K + K~), and 77 — > K + K~7t + tt~ that have no K° s candidates, two different 
sets of selections are applied to the Kg candidates with high (low) momentum, 
i.e. with momentum > 1.5 GeV/c (0.5-1.5 GeV/c): dr is required to be larger 
than 0.02 (0.03) cm for both charged tracks; the 7T + n~ vertex is required to 
be displaced from the IP by a minimum transverse distance of 0.22 (0.08) cm. 
The mismatch in the z direction at the Kg vertex point for the 7r + 7r~ tracks 
must be less than 2.4 (1.8) cm; the direction of the pion-pair momentum 
must also agree with the direction from the IP to the vertex to within 0.03 
(0.1) rad. To evaluate the background and calculate efficiencies, we use a 
Monte Carlo simulation (MC) of the detector response based on GEANT3 
[9]. The TREPS code [10] is used for 77 — > K S K S event generation and the 
background 77 — > 2(7r + 7r~), 77 — > 2(K + K~) event generation. From MC 
simulation, with the described K s selection above the K s signal efficiency can 
reach ~ 80% while the background is reduced by a factor of 10 5 . Thus the four- 
track backgrounds can be eliminated efficiently after our event selection. The 
resolution in the reconstructed Kg mass is 4 MeV/c 2 , and only candidates for 
which |M(7r + 7r~) — m K o \ < 13 MeV/c 2 are selected. Finally, we require that 
the sum of the transverse momentum vectors of all tracks in the cm. frame 
of the e + e~ beams, |Spf e |, be smaller than 0.1 GeV/c (momentum balance). 
W is calculated from the invariant mass of the KgKg pair, and | cos 6* | is 
obtained from the Kg scattering angle with respect to the incident axis of 
the electron in the 77 cm. frame, which approximates the direction of the 
incoming photon. Figure 1 shows the tt + it~ invariant mass spectra after stage 
I and stage II selection. After applying the above selections, we find 981 KgKg 
candidates in the range 2.4 GeV< W < 4.0 GeV and | cos6>*| < 0.6. The W 
distribution is shown in Fig. 2. Clear signals for the Xco an d Xci resonances 
are observed. 



3 Background subtraction 



The background contamination from events where additional particles accom- 
pany the two detected Kg mesons - so-called non-exclusive backgrounds - 
should be also estimated. Because of the available phase space, such events 
are expected to have a | Sp| e | distribution that is close to zero at |Sp^ e |=0 
and increases with |Sp^ e |. This feature is verified in the 77 — > KgKgir (which 
is the dominant background) MC and data sample, where the MC sample is 
generated by using GGLU code [11]. We assume that the |Spf e | distribution 
of the non-exclusive background can be parameterized by 
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0.48 0.49 0.5 0.51 

M(tcV) (GeV/c 2 ) 



0.52 



Fig. 1. The 7r + 7r invariant mass spectrum for Kg candidates after stage I (dotted 
histogram) and stage II (points with error bars) K$ selection. Here events are se- 
lected in the range W=2.3-4.5 GeV and |Spf| < 0.25 GeV/c, where W and |Spf| 
are calculated by assuming all tracks are charged pions. 
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Fig. 2. KgK® mass spectrum before background subtraction. The curves show the 
fit result described later in the section on Xcj resonances. 




cx, x < 0.05 (GeV/c) 
ax + b, x > 0.05 (GeV/c) 



constrained by 0.05c = 0.05a + 6. We fit the function f(x) to the difference 
between data and signal MC distributions which is normalized to the data be- 
low 0.03 GeV/c where the background contribution is negligibly small (Fig. 3). 
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Using the data sample with |Spf e | =0.5-1.0 GeV/c we verify that there is no 
9* dependence of the shape. Using the fit results, the estimated background, 
which is 4.1 ± 0.1%, 3.6 ± 0.2%, and 2.6 ± 0.3% for 1^=2.4-2.6 GeV, 2.6- 
2.8 GeV, 2.8-3.3 GeV, respectively, is subtracted in each W bin (the errors 
are statistical only). For W=3. 6-4.0 GeV, the background is set to zero since 
the data sample is too small to apply the procedures described above. 

Finally, 952 signal events remain in the signal region | Sp| e | < 0.1 GeV/c after 
background subtraction. 




°0 0.05 0.1 0.15 0.2 

EpP^GeV/c) 

Fig. 3. |Xpf e | distribution for K^Kg candidates. The dotted histogram and points 
with error bars indicate the distribution of events before and after background 
subtraction, respectively. The dot-dashed line is the background distribution, which 
is obtained from the fit to the difference between MC and data. The solid curve 
shows the signal MC distribution, which is normalized to the number of signal 
candidates in the three leftmost bins. The arrow indicates the upper boundary of 
the |Spf e | requirement for the signal. 



4 Cross sections of the process 77 — > K^Kg for 2.4 GeV < W < 4.0 GeV 



The differential cross section for two-photon production of the final state X 
in electron-positron collisions is given by 



■{W, \cos6 |; 77 ->X) = /ttt-\ a tt t a i ' ' i 7777 i n*\\ n @ ) 



d|cos0T " ' ' C rr (W)AWA\co8e*\e(W,\cose*\)C i 

where and e denote the number of signal events after background sub- 
traction and the product of detection and trigger efficiencies, respectively. The 
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integrated luminosity of this experiment, £j nt , is 397.6 fb 1 and is determined 



with a systematic uncertainty of 1.4 
function of W, is defined by 



/o. 



The luminosity function £ 77 , 



e ' e 



-x) 



(4) 



The efficiencies e(W, | cos#*|) are obtained from MC using the TREPS code 
[10] for 77 — > KgKg event generation. The TREPS code is also used for the 
luminosity function determination. Trigger efficiencies are determined from the 
trigger simulator. The typical values of the detection and trigger efficiency are 
5-19% and 90-95%, respectively, and grow with increasing W and decreasing 
| cos 6* | . Differential cross sections normalized to the cross section integrated 
over the range | cos#*| < 0.6 (o"o) in different W bins are shown in Fig. 4(a). 
The angular distributions are consistent with both BC and DKV predictions 
up to | cos9*\ = 0.5. 

The angular distributions, a^da/d] cos(9*) |, in the Xco and \ci regions — 
m(Xco)| < 66 MeV/c 2 , \W - m(x<a)\ < 36 MeV/c 2 ) shown in Figs. 4(b,c) 
are in good agreement with those expected for the decays of the spin zero 
and two particles. The total cross section o"o as a function of W is shown in 
Fig. 5(a) and listed in Table 1. The values of the total cross section for the 
range, W / =3.3-3.6 GeV, where the contribution from charmonium states is 
large, are omitted. 



W(GeV) 


iVev 


o"o, nb 


2.4-2.5 


226.3 ± 15.4 


0.0816 ±0.0056 ±0.0070 


2.5-2.6 


195.6 ± 14.3 


0.0671 ±0.0049 ±0.0057 


2.6-2.7 


137.9 ± 12.0 


0.0488 ± 0.0042 ± 0.0042 


2.7-2.8 


81.9 ±9.2 


0.0307 ±0.0034 ±0.0027 


2.8-2.9 


46.8 ±6.9 


0.0178 ± 0.00261^021 


2.9-3.0 


31.1 ±5.7 


0.0131 ± 0.00241Hm6 


3.0-3.1 


21.4 ±4.7 


0.0084 ± 0.0018±8;88io 


3.1-3.2 


10.7 ±3.3 


0.0046 ± 0.00141™ 


3.2-3.3 


10.7 ±3.3 


0.0039 ± 0.0012±8;ggg| 


3.6-4.0 


5.0 ±2.2 


0.0006 ± 0.0003iHooe (< 0-0013 at 90% CL) 



Table 1. Signal yields (iV cv ) and total cross sections (co) for the process 
77 — > KgKg in the angular range |cos#*| < 0.6. The first and second errors 
are statistical and systematic, respectively. 
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0.2 0.4 

|COS9*| 



Fig. 4. (al-a4) The angular distribution of the cross section, &Q l da/d\ cos 9*\, in 
different W ranges. The solid curves are 1.227 sin" 4 9*, which is the prediction of 
DKV. The dotted curves are the prediction of BC. (b) The angular distribution in 
the XcO region; the dotted curve shows a flat distribution (J=0); (c) The angular 
distribution in the Xc2 region; the dotted curve shows the helicity 2 distribution 
(oc sin 4 9*). The errors indicated by short ticks are statistical only. 



5 Systematic errors 



The dominant systematic errors are summarized in Table 5. We assign 4% 
to the uncertainty from trigger, which is determined by comparing the trig- 
ger efficiencies in the data sample and trigger simulation. The uncertainty 
of K s reconstruction efficiency is estimated by comparing the ratio of the 
number of 77 — > K S K S events with both K s mesons satisfying the selection 
requirements and that with only one K s satisfying the requirements in data 
and MC samples. We take the efficiency difference between the data and MC 
77 — > K S K S sample, which is 4.4% for one K s . The uncertainties in the 
background subtraction are estimated by fitting the background shape in the 
I Sp| e I distributions using second-order polynomial functions and comparing 
the background fractions obtained to those described above. The differences 
between the two calculations are taken as the corresponding systematic error 
in each energy range and are 2.0%, 2.0%, and i|i% for W=2A-2.Q GeV, 2.6- 
2.8 GeV, 2.8-3.3 GeV, respectively. For W=3. 6-4.0 GeV, we conservatively 
assign the number of observed events as the systematic error in the back- 
ground. The 3.4-5.0% systematic error for the luminosity function in the range 
W=2. 4-4.0 GeV in Ref. [10] is determined from comparison of the kinematic 
distributions for the two-photon system in events generated with TREPS to 
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those from a QED calculation that includes all order a 4 diagrams [12]. The 
total W^-dependent systematic error is (8.5-12.1)%. 



Source 



Error, 



Trigger efficiency 4 

Luminosity function 3.4-5.0 

Background (for non-resonant analysis) 2.0-8.4 
K s reconstruction (per Kg) 4.4 
Integrated luminosity 1.4 



Total 



8.5-12.1 



Table 2. Summary of systematic errors 



397.6 fb" 1 




0.15 



0.1 



0.05 



1 a (K^)/a (K + 
I f (IcosB l<0.6) 

tjj>. 



(b) 



DKYSUCMJimit... 



in" 4^ 

BC.predidion. 



2.5 3 3.5 4 2.5 3 3.5 4 

W(GeV) W(GeV) 



Fig. 5. (a) Total cross sections for 77 — > K S K S in the cm. angular re- 
gion I cos 9*\ < 0.6. Here n is the IT-dependence (oc W~ n ). (b) The ratio 
a (K s K s )/a (K + K~) versus W in | cos 9*\ < 0.6, where the K + K~ data are taken 
from the Belle measurement [4]. The dotted line is the DKV prediction with the 
flavor symmetry assumption; the dashed and dashed-dotted lines are the BL and 
BC predictions, respectively. The two sets of error bars show the statistical and 
combined statistical+systematic errors, respectively. 



6 Discussion 



The leading term in QCD calculations [1,2] predicts a ~ W~ G dependence of 
the cross sections cicx/ci cos #* (77 — > MM). However, the fit to the data in 
the range W=2A-A.O GeV gives a W-dependence (<r oc W~ n ) of n = 10.5 ± 
0.6 ± 0.5, where the first error is statistical and the second is systematic. We 
conservatively estimate the systematic error on n by artificially deforming the 
measured cross section values assuming that the systematic errors are strongly 
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correlated point-to-point, as in Ref. [4]: we shift the <7o values at the two end 
bins by ± 1.5 and =(=1.5 times the systematic error, respectively, whereas each 
intermediate point is moved so that its shift follows a linear function of W 
times its systematic error. The average of the observed deviations in n from 
its original value is taken as a final systematic error. The value of n indicates 
that, unlike 77 — > -k + -k~ and 77 — > K + K~ [4], the current values of W are 
not yet large enough to neglect power corrections in 77 — > K^R®, which are 
not taken into account in the BL and BC predictions. 

The ratio ao(KgKg)/ao{K + K^) shown in Fig. 5(b) decreases from ~0.13 
to ~0.01 with increasing W. This energy dependence is inconsistent with the 
DKV prediction that the ratio should be ~ 2/25 in the SU(3) symmetry limit. 
Furthermore, it is difficult to explain the experimental result with the handbag 
model even if the effect of SU(3)-symmetry breaking is taken into account 
[6,13]. This indicates that the handbag model needs significant corrections. 

Since the experimental values of the ratio ao(KgKg)/ao(K + K~) approach 
the BL and BC predictions at the highest measured energies W ~ 4 GeV, the 
leading term QCD calculations [1,2] may become applicable for a(KgKg) at 
not much larger values of W. 



7 The two-photon decay width of Xcj resonances 

Measurements of 77 — > K®Kg can also provide more precise results [14] for the 
two-photon decay widths and branching fractions of the charmonium states 
since the continuum background is strongly suppressed. By fitting the contin- 
uum M(KgKg) distribution to an exponential distribution and parameterizing 
the charmonium peaks with a Breit-Wigner function for the Xco and Gaussian 
function for the narrow \c2 with the masses and widths floating, 134 ± 12 Xco 
and 38 ± 7 Xc2 events are observed. The masses and widths obtained from the 
fit taking into account the detector resolution are consistent with the PDG 
values. The Xco{Xc2) statistical significance is 22.7<r (11.2<r), where a is a stan- 
dard deviation. The statistical significance of the signals is obtained from the 



yj — 2 ln(Lo/-L m ax) values, where Lo(max) is the likelihood without (with) the 
signal contribution, with the joint estimation of the three parameters (mass, 
width, and yield are determined simultaneously). The two- photon decay width 
of the Xco or Xc2 can be obtained using the formula 




r 77 (x c j) x B(xcJ 




Ym 2 



4(2J + l)ir 2 C^(m)eB 2 (K° s 



7T + 7T )£ 



hit 
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where Y and m are the yield and mass of the Xcj charmonium state, re- 
spectively. The quantity e denotes the product of the detector efficiency, trig- 
ger efficiency, and angular acceptance for the resonant decays. In addition 
to the sources of systematic errors listed in Table 5, the errors in the yield 
are 2.3% and 2.4% for the Xco and Xc2, respectively. For Xc2 events we as- 
sume a pure helicity 2 state in MC generation following the previous measure- 
ment [15] and theoretical expectations [16,17]. The directly measured values of 
the product T^( X cj)B(Xcj -> K° S K° S ) are 7.00 ± 0.65 ± 0.6 eV for the Xco and 
0.31 ± 0.05 ± 0.03 eV for the Xc2- Using the results of our previous measure- 
ment of K + K~ and 7r + 7r~ production in 77 collisions [4], we determine the 
ratios B(K° S K° S )/B(K + K-) and B{K s K q s )/B{ti + ^) for the X co and X c2, in 
which some common systematic errors cancel. Here B(tt + tt^), B(K + K~) and 
B(KgKg) are the branching fractions for the Xcj decay to the corresponding 
final state. Using the world-average values of the branching fractions B(xco ~~ > 
K° S K° S ) = (2.8 ± 0.7) • 10~ 3 and B( X c2 -> K° S K° S ) = (6.7 ± 1.1) • 10~ 4 [14], 
from the products of the widths and branching fractions given above we can 
extract the values of the two-photon width that are shown in Table 7. The 
notation for. indicates the systematic uncertainty from the branching fraction 
°f Xcj — > Ks-Ks- ^ c an be seen that for both the Xco and Xc2 the value of 
B(KgKg)/B(K + K~) is compatible with 0.5 as expected from isospin symme- 
try. The values of the two-photon widths of the Xc0(2) charmonia are consistent 
with those obtained from their total width and the branching fractions for de- 
cay to two photons in Ref. [14]. 



Resonance Xco Xc2 

T. n B(K s K° s ), eV 7.00 ± 0.65 ± 0.71 0.31 ± 0.05 ± 0.03 

B(K° S K° S ) /B(K+K~) 0.49 ± 0.07 ± 0.08 0.70 ± 0.21 ± 0.12 

B(K s K° s )/B(7T + 7i-) 0.46 ±0.08 ±0.07 0.40 ± 0.10 ± 0.06 

T 77 , keV 2.50 ± 0.23 ± 0.23 ± 0.62 (6r.) 0.46 ± 0.08 ± 0.04 ± 0.08(6r.) 

Table 3. The products of the two-photon width and the branching fraction, ratios of the 
branching fractions, and two-photon widths for the XcO an d Xc2- The notation br. indicates 
the systematic uncertainty from the branching fraction of Xcj KgKg 



8 Conclusion 



Using a 397.6 fb -1 data sample accumulated with the Belle detector at KEKB, 
the cross sections of the process 77 — > K° S K^ have been measured for the 
first time in the W range from 2.4 GeV to 4.0 GeV with |cos#*| < 0.6. 
The overall IV-dependent systematic uncertainty is 8.5-12.1%. The measured 
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W- dependence (er oc W~ n ) of 77 — > K° S K° S is n = 10.5 ± 0.6 ± 0.5 from 
a fit to the data with W=2. 4-4.0 GeV, indicating that, unlike 77 — > 7r + 7T~ 
and 77 — > K + K~, the W values up to 3.3 GeV are not sufficiently large to 
apply the leading term BL and BC predictions to 77 — > K^Kg. The angular 
distribution in the range | cos#*| < 0.5 is consistent with both BC and DKV. 
The ratio 0"o(77 — > K^Kg) / '00(77 ~^ K + K~) decreases rapidly from ~0.13 to 
~0.01 with increasing W in contrast to the expectation from the DKV model. 
Since the measured values of the cross section ratio approach the BL and BC 
predictions in the highest energy bin, 3.6-4.0 GeV, this may indicate that the 
leading term QCD calculations for (7(77 — > KgKg) are already applicable at 
W values larger than ~ 4 GeV. In addition, the products of the two-photon 
decay width and branching ratio to Kg Kg for the Xco an d Xc2 are found to be 
7.00 ± 0.65 ± 0.71 eV and 0.31 ± 0.05 ± 0.03 eV. 
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